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SUMMARY 


A theoretical analysis is presented for the harmonic noise of high speed, open 
rotors. Far field acoustic radiation equations based on the Ffowcs-Williams/Hawkings 
theory are derived for a static rotor with thin blades and zero lift. Near the plane of 
rotation, the dominant sources are the volume displacement and the pu^ quadrupole, 
where u is the disturbance velocity component in the direction of blade motion. These 
sources are compared in both the time domain and the frequency domain using two- 
dimensional airfoil theories valid in the subsonic, transonic, and supersonic speed 
ranges. For nonlifting parabolic arc blades, the two sources are equally important 
at speeds between the section critical Mach number and a Mach number of one. How- 
ever, for moderately subsonic or fully supersonic flow over thin blade sections, the 
quadrupole term is negligible. It is therefore concluded for thin blades that significant 
quadrupole noise radiation is strictly a transonic phenomenon and that it can be sup- 
pressed with blade sweep. Noise calculations are presented for two rotors, one simu- 
lating a helicopter main rotor and the other a model propeller tested at United 
Technologies Corporation. For the latter, agreement with test data was substantially 
improved by including the quadrupole source term. 

INTRODUCTION 

In the 1970's there has been a renewed interest in noise of open rotors, not only 
for helicopter application, but because of a development program in the United States 
for fuel conservative propulsion systems, Ffowcs -Williams and Hawkings provided a 
theoretical basis (ref. 1) for analyzing this type of problem by showing that, in prin- 
ciple, the noise can be calculated exactly if certain aerodynamic quantities (source 


This paper was originally presented at the Spring Meeting of the Institute of Acoustics, 
Cambridge University, Cambridge, England, April 7, 1978. The manuscript has been 
submitted for possible publication in the Journal of Sound and Vibration and is repro- 
duced herein by permission of the editor of that journal. 
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terms) are known on and around the blades. The surface sources are the volume dis- 
placement and the blade surface forces. The volume sources are quadrupole terms 
representing shear stress in the air. The propeller noise problem is an ideal applica- 
tion of the theory because the blades are thin and, in forward flight, the sources are 
essentially steady in the blade coordinate system. 

Because the surface source terms are linear, considerable progress has been 
made in developing techniques for computing their noise (refs. 2, 3, 4). However, the 
quadrupoles have been neglected because they are second order in the disturbance 
velocities, which makes them appear small. This appeared justified because the linear 
theories worked reasonably well for moderate speed propellers. Also, sonic boom 
theory, which uses only the linear volume displacement and loading sources, is suc- 
cessful at supersonic aircraft speeds. However, Kitaplioglu and George (ref. 5) re- 
marked in a recent paper that the linear theories consistently underpredict noise from 
rotors operatic^ at the transonic speeds which are currently of interest. An example 
of this is the transonic propeller in figure 1 which was run at United Technologies 
Research Center. Figure 2 shows test versus linear theory from ref. 4. Maximum 
noise along the fuselage occurs near the plane of rotation, where noise levels predicted 
from the surface sources are as much as 5 dB lower than the data. After considering 
other mechanisms, the quadrupole source was investigated as an explanation for this 
underprediction. 

The purpose of this report is to illustrate the role of the quadrupole sources in 
noise radiation from open rotors using the simplest possible meaningful example. To 
accomplish this, far field radiation equations are derived for a static, nonlifting 
rotor. The equations are cast into a form which permits direct comparison of the vol- 
ume displacement and quadrupole sources inside the radiation integrals. The source 
terms are then evaluated using aerodynamic techniques which are valid through the 
subsonic, transonic, and supersonic ranges. Finally, some sample noise calculations 
are presented. 


SYMBOLS 

b airfoil chord 

B number of blades 

Bjq = b/D chord-to-diameter ratio 
Co ambient speed of sound 

D propeller diameter 

fj i^ component of force/unit area exerted by airfoil on fluid 

G Green’s function (eq. 2) 

h thickness distribution (figure 3) 


340 



H 

1» j 

Jn 

k 

m 

Mj, 

n 

P 

Pn 

PVn 

PQn 

r 

To 

Trp 

R 

s 

t 

T 

T • ■ 
ij 

u 

U 

V 

Vn 

X 

X 

y 

yi 


hA>» normalized thickness distribution (figure 4) 

1, 2, 3 Cartesian coordinate indices 
Bessel function 

specific heat ratio, 1.4 for air 
harmonic of blade passing frequency 

S2 — , section relative Mach number 

Co 

n — , tip rotational Mach number 
Co 

mB, harmonic of shaft frequency 

acoustic pressure 

complex Fourier coefficient of p 

volume displacement portion of Pjj 

quadrupole portion of 

distance from origin to observer point 

distance from origin to source point on blade 

propeller tip radius, D/2 

I iT - ^ I , distance from source point to observer 
source surface in eq. 1 
observer time 

integration limit for source time (r) integration 

pujuj, quadrupole (shear stress) source 

component in chordwise direction of disturbance velocity 

component in i^ direction of disturbance velocity 

nro , local blade section speed 

component normal to chord of disturbance velocity 

normal component of airfoil surface velocity 

(x, y, 0), observer coordinates (figure 3) 

y/b, normalized chordwise coordinate 

observer distance from propeller axis (figure 3) 

source coordinates yj, y 2 , ys 
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l/b, normalized coordinate perpendicular to chord and to radius 

— , normalized radial coordinate 

source coordinate normal to chord (figure 3) 
source coordinate in chordwise direction (figure 3) 

Dirac delta (impulse) function 

angle from propeller axis to observer point 

source volume (volume exterior to blades) 

density 

ambient density 
source time variable 

integrated volume displacement source strength in frequency domain 
(eqs. 31, 34) 

integrated quadrupole strength in frequency domain (eq. 33) 
integrated strength of pu^ quadrupole in frequency domain 
radian frequency of sound 
2 n times shaft rotation frequency 


I 


indicates vector, as in x 

indicates differentiation with respect to argument, as in h^ 


ACOUSTIC THEORY 

The starting point for the analysis is Goldstein's version of the acoustic analogy 
(eq. 3.6, ref. 6) from which the Ffowcs -Williams/Hawkings formulas (ref. 1) can be 
derived: 


p’(x,t) = 


T 

- f f 

Co J-T Js{T) 


3G ^ dG 


dS(y) dr 


+ 




( 7 -^ oyiayj 


a2G ^ 

7 dy dr 


( 1 ) 


This equation gives the disturbance density p' exactly for known values of the 
source terms. The sources to be evaluated on the moving surface S (t) are the nor- 
mal surface velocity Vn (taken to be positive outward, the opposite of Goldstein's 
convention) and the surface force components fj. The quadrupole source is the 
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Lighthill stress tensor Tjj which is to be evaluated in the volume around the blades 
v{t). The source time (t) integration is over a range - T <r<T large enough to 
include all signals from the source region which arrive at the observation point at 
time t. G is the Green's function 


6(t-r-R/co) 

47tR 


( 2 ) 


where R = [jT-yj is the distance between source point ~y and observer point jT. 

For a thin, nonlifting propeller blade, the surface forces fi can be neglected. 
Also, we make the usual approximations that the acoustic pressure 


p = c2 p' (3) 

where c^, is the ambient sound speed and 

Tij = pu^Uj ^4^ 

where p is the density and ui is the disturbance velocity component in the i^^ direc- 
tion. The pressure term in Tjj is neglected for now but can be added at any time in 
the analysis. For the static propeller, we approximate the Cartesian source coordi- 
nates with locally orthogonal curvilinear coordinates as shown in figxure 3: 


Yl = 7 Y2 = ^ ^3 = 


The observer coordinates are x, y, 0. 


The distance R is given by 

R = \j {'x. + - 2yro cos i—\ 

' ^ O ' 

In the far field Greenes function 


( 6 ) 


Q ^ 6(t-r-R/co) 
47T r 


(7) 


only the terms of order |/r and ro/r in a series expansion of R are retained, giving 


R 


r + ^ cos d - 


sin B cos 



where B is the angle of the observer from the propeller axis. 


(8) 
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Since the blades of interest are, thin, we consider the surface sources to act on the 
mean chord line ^ = 0 and neglect the volume occupied by the blade. This is easily 
justified, for example, with blades 2% thick and a 0. 15 chord-to-diameter ratio opera- 
ting at a tip Mach number of 1. 0 because the maximum error in point of action of the 
sources is less than the wavelength of the 2000 ^ harmonic of shaft rotation frequency. 

To evaluate the surface integral at r = 0, we note that Vn dS = — U h' ( 7 ) dy dro 

for each surface of the airfoil. Here U = n ro is the local blade section speed, h(y) 
is the airfoil thickness distribution shown in figure 3, and the prime on h denotes dif- 
ferentiation with respect to the argument. The only time dependence of the sources 
considered here is that due to convection at speed U in the negative 7 direction, so 
the source behavior is given simply by h (y+ Ur) and Tij (7+ Ur). 


From these arguments, eq. 1 reduces to 

— /■“ sr 

P(x,t) = / / / -P„Uh ' (y + Ur) dy dr„ dr 

J-CO Jq J-TTV^ 


^00 /.oo r 2 

l-JJ-nr.L Irfr * 


( 9 ) 


d^ dydr. dr 


where the integrations over upper and lower airfoil surfaces are lumped together using 
y as an integration parameter and rrj, is the blade tip radius. T has been set to 00 , 
assuming the t integrals will converge. The sum on i and j has been performed 
neglecting any radial velocity disturbances so that the remaining quadrupole sources 
are 


Til 

(10) 

Ti2 

(11) 

^22 

(12) 


where u and v are as shown in figure 3. The great simplification brought to eq. 9 
by the thin blade approximation is that the y and | integration limits no longer 
depend on r . Thus, the integrations may be performed in any order. Equation 9 
gives the noise caused by one blade; the other blades are accounted for by superposi- 
tion. 


The-derivatives of the Green's function with respect to r and 7 are shifted 
temporarily onto the source functions using integration by parts. 
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p{x,t)=J I I Po h" (7+ UT) G dr dr„ dr 

J-K Jo J-TTTo 


/*" r°° /•“ r " ' 

+ / / / / (Y+ Ur) G - 2 Tj2 (r+ UT) 

J-qoJq L 


a| 


+ T22 (Y+ UT) 


8^g ' 

a|2. 


dy dro dr 


(13) 


Because h’ is discontinuous at the blade leading and trailing edges and because 
the Tjj's are discontinuous across shocks, their derivatives must, at this point, be 
considered as generalized functions in the sense described by Farassat <ref. 7). The 

9 G 

differentiations will shortly be removed from the sources. The derivatives — and 

q2q 

can be evaluated explicitly from eqs. 7 and 8: 


aG 

a? 


cos & 

47TCor 


5' (t- t-R/cJ 


(14) 


a^G 

a 


cos ^6 
4 7T c r 


6" (t-r-R/cJ 


(15) 


Substitution of these into eq. 13 eliminates derivatives with respect to | . The 
derivatives with respect to the argument of the delta function can be removed using in- 
tegration by parts on the r integral. 


p(x, t) 


] f f " r 

< nr I I f 

1 /•“ r " u CO! 

/ / / / T,, (y* Ut) + 2 : 

47rr/// / 

•^-““'0 •'-nr^J-oc L 


PgU h (ytUr)(5(t-T - R/c„) dv dr„ dr 


T,2 0'+ Ut) * 


cos^ 0 


T 22 (v* UT) 


(16) 


5(t-T-R/c„) d| dy dr„ dT 


Time Domain Radiation Equations 


The T integration is now trivial, simply replacing the arguments 7+ Ur by 
y+ Ut - MrR where Mr = U/c = Q.v^/Co is the section relative Mach number. 


p(x, t) = 


1 o 

4in / / Po h” (y+ ut - Mj.R) dy dr„ 

•X J- nr„ 


* J J [^’’^11 2 T^2(y+ Ut-MpR) + cos^e T 22 (y+ Ut-Mj,R)j d^dydr^ 


(17) 


The derivatives are removed from the source functions by noting that 

h (y+ Ut - MrR) = U^h (7+ Ut - MrR) and that can be moved outside the 

at'^ at2 

integrals because the limits are independent of t . 
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= 47F7 P.,h(Y+Ut-Mj,R) dv dr„ 

f T22 ( 7 +m-MrR) d? dY dr, 

dtryQy__j, j_^ u o of 


(18) 


The y integrations in eqs. 17 and 18 run from -TTr^ to TTr^. However, 
h (y+ Ut - Mj.R) is zero except for values of y satisfying 

2 TTNro - -| < y + Ut - M^R < 2 TrNro + ^ 


(19) 


where N is an integer. Finding these values of y is equivalent to finding the acous- 
tic planform or retarded blade location as described in ref. 4. 

The first integral in eq. 18 is the traditional thickness noise which is the direct 
result of volume displacement at the surface of the blade. This surface integral over 
the rotor disc was derived previously by Hanson (ref. 4) and was shown to give the 
same result as Farassat's theory (ref. 2). The second integral in eq. 18 is the quad- 
rupole thickness noise which is also a result of volume displacement but is accounted 
for by integrating the associated stress terms T^^* T^i 2 * '^22 the volume 

surrounding the rotor disc. The quadrupole term has not been calculated correctly in 
previous work. 


Equation 18 has been presented in a form which shows that the T^2 quadrupole 
(integrated over | ) radiates exactly like the volume displacement source with regard 
to frequency and directivity dependence. Since T 12 und T 22 ure multiplied by cos 0, 
we need only compare chordwise distributions of 


h 

b 


versus 



( 20 ) 


to compare contributions of the two sources to noise near the plane of rotation. This 
is done in a later section entitled "Aerodynamic Evaluation of Quadrupole Source 
Term. " 

The time domain equations can be transformed to blade fixed coordinates by sub- 
stituting 71 = y + Ut - Mj.R. Since this might be an advantage for subsonic tip speed 
rotors, an example is given in the appendix. 

Frequency Domain Radiation Equations 

It is also useful to compare these sources in the frequency domain, which is done 
by calculating complex Fourier coefficients from either eq. 17 or eq. 18 according to 
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>n <^) = 

*'0 


2n/Q 


p (jT, t) e dt 


( 21 ) 


The volume displacement term from eq. 18 becomes 

2 TC 

,Trj. _7rr 


Pvn( 5 J‘) = - f f ° ^ " h (V+ Ut - M^R) 6^“"* dt dV dr, (22) 

Jo J-TTVo Jo 


where the factor in the frequency domain comes from — ^ in the time domain. 

at2 

We change variables with 


Xb = T + Ut-MrR 


(23) 


Then, because the surface source is assumed to act on the mean chordline ^ = 0, 
eq. 8 gives 


Xb y r ro . ^ 
t = 7 T — - — — + — - — sin e cos 

f^rQ ^ TC Q Cq Co 


e) 


This results in 


(24) 


Pvn(^) = - 


Pc,n 2 fi 2 Co 

e 

4 nr 


nr r 


: rrp ^ 1 

X ii: 


in- 


bX 


h(Xb) e ° dX 




V S2r v 

2nro -in— - in sin 6 cos — 

e dr 


dr. 


(25) 


The y integration gives Bessel functions so that eq. 25 reduces to 


Pvn(x) = - 


2^2 in 
Pon n 

47rr 




nb 

i— X 

h (bX) e dX dro (26) 


which is the same result derived by Hawkings and Lowson (ref. 3). Equation 26 Is now 
normalized using 


z = — 

rT 


(27) 
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(28) 
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nrrp 


Co 
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D 



(29) 


and the thickness function H (X) = h (bX)/b shown In figure 4 to arrive at the final form 
for the volume displacement noise component; 



where 





2nB 


D 


H(X) e 


dX 


(31) 


is the Fourier transform of the airfoil thickness distribution. The same manipulations 
for the quadrupole term in eq. 18 give 


^Qn(^) - ~ Po Co 


2 




Jjj (nZ M-p sin 9) 


X + 2 Z cos e 


tjfl2 + cos^e J/^22) 


(32) 


where Y = l/b and 



y 2nBj) 


(”?)■ 

r y 21nBD Mt cos 9 Y , ‘ Z ^dx 

(33) 

\ z / 



is the integrated quadrupole source term. The chordwise integration is over a range 
± Xj equivalent to ±$2ro. Because eqs. 30 and 32 are in the same form, the two 
sources can be compared in the frequency domain by comparing the source terms 
xlfy and ij/ii . This and the time domain comparison are the subjects of the following 
section. 
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AERODYNAMIC EVALUATION OF QUADRUPOLE SOURCE TERM 


It is well known from experimental studies that noise of high speed rotors is max- 
imum near the plane of rotation. Since cos 0 = 0 in the plane of rotation, the pre- 
ceding acoustic analysis shows that the volume displacement and quadrupole thickness 
sources can be compared by calculating 


and 





2mBBr) 

i X 


H(X) e 


dX 


(34) 


>Uii 



2mBBD 

e 


X 


dY dX 


(35) 


where n has been replaced by mB and m is the harmonic of blade passing fre- 
quency. The streamwise velocity perturbation is u and p is the local density. 
The chordwire integration is shown with an infinite range because the source is 
now considered to be "unwrapped" to allow use of two dimensional airfoil theory for 
evaluation by These source terms have been evaluated for biconvex para- 

bolic airfoils. 


H(X) = (1 - 4 X^) (36) 

b 

and the results are presented below. 

Aerodynamic Theory 

The airfoil flow field was studied using different computation schemes for fully 
subsonic, transonic, and fully supersonic flow as described herein. All calculations 
were two dimensional and assumed small disturbances so that the density ratio could 
be calculated from the small perturbation approximation for isentropic compressible 
flow 


_P 

Po 




"(U + u)2 + v^ 


U^ 



1 - Mi 


u 

TT 


where k is the ratio of specific heats. 


(37) 
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For fully subsonic flow, velocity distributions were calculated using the Spreiter 
and Alksne local linearization method as given in ref. 8. 

For fully supersonic flow, e3q>licit first order solutions for the supersonic flow 
field of a parabolic arc airfoil were taken from Caughey (ref. 9). The calculated shock 
wave positions for a 2% thickness ratio airfoil at a Mach number of 1. 15 are plotted in 
figure 5. Because of the spreading shock waves, the quadrupole noise originates from 
an increasingly noncompact region at increasing heights above the airfoil. Streamwise 

perturbation velocities near the surface are predicted to vary inversely as J - 1 , 

as is also predicted by linearized supersonic theory. Thus, the ratio of quadrupole to 
volume displacement acoustic pressure is predicted to increase as supersonic Mach 
number is decreased and the largest value of this ratio should occur in the transonic 
regime. 

For transonic flow, elaborate digital computer programs are available for detailed 
calculation of transonic flow past two-dimensional airfoils. Rather than modify one of 
those programs to calculate the spatial distribution of perturbation velocity, an older 
approximate method for transonic flow was modified by use of information available 
from the newer, more rigorous solutions. This modified solution is also based on the 
local linearization method developed by Spreiter and Alksne (ref. 8). In that method, 
the streamwise velocity gradient in a locally subsonic portion of the flow field is ob- 
tained from its value for Incompressible flow by 

Here, the quantify in square brackets is one minus the local axial component of Mach 
number squared. Similarly, streamwise velocity gradient in a locally supersonic flow 
field is obtained from its value as calculated from linearized supersonic theory at a 
free stream Mach number yj~2 by 

1 

A. (v) = [m? - 1 - (k+l) M2 ^ 

Also, streamwise velocity gradient in an accelerating near-sonic flow was shown to be 
equal to its value at the same location at a free stream Mach number of one, divided 
by free stream Mach number to the 2/3 power. 

Equations were given by Spreiter and Alksne in ref. 8 for calculating the stream- 
wise perturbation velocity within these three regions, assuming that the value of this 
velocity was known at one chordwise position. This velocity, and therefore the con- 
stant of Integration when integrating the above velocity derivatives, was known only for 
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the cases of fully subsonic, near sonic, and fully supersonic flow. Fortunately, numer- 
ical solutions for the flow field of a nonlifting parabolic arc airfoil are available for 
transonic Mach numbers less than one (Murman and Cole, ref. 10) and greater than one 
(Murman, ref. 11). These numerical results can be utilized to evaluate the variations 
of any flow property of interest throughout the transonic range. In particular, it had 
been noted in ref. 8 that the equation for perturbation velocity in accelerating near- 
sonic flow had a denominator that went to zero at the sonic point on the airfoil surface. 

To obtain valid results, the numerator also had to go to zero such that their ratio re- 
mained locally constant. The anal 3 dical solutions for fully subsonic flow, fully super- 
sonic flow, and free stream Mach numbers very near one and the numerical results in 
refs. 10 and 11 for the intermediate regions were utilized to prescribe the streamwlse 
perturbation velocity at this chordwise position (for a parabolic arc airfoil, the quarter- 
chord). 

Velocity gradients calculated from eqs. 38 and 39 should be reasonably accurate if 
the local Mach number is not near one. As local Mach number approaches one, the cal- 
culated gradients approach infinity. Therefore if one starts with a prescribed locally sub- 
sonic velocity at some chordwise location and numerically integrates eq. 41 in the down- 
stream (increasing local velocity) direction, incorrect results would be obtained at near - 
ly sonic local Mach numbers. As with the calculation method described by Fink (ref. 12) 
for axisymmetric transonic flow, the numerical integration has been modified to use 
velocity gradients calculated for accelerating near-sonic flow at the chordwise position 
where it gives a smaller gradient. Further downstream where the flow was locally 
supersonic, the velocity gradients calculated for accelerating transonic flow and for 
locally supersonic flow were compared, and the smaller value was used in numerical 
integration. 

Shock wave position was prescribed by use of a curve fitted to the positions given 
in ref. 10 for this airfoil section. Shock wave strength in the presence of a boundary 
layer was obtained by an approximation to the data correlation given by Sinnott in ref. 

13. The variation of perturbation velocity with distance above and below the airfoil was 
calculated from the Spreiter and Alksne theory (eqs. 63 and 64 of ref. 14) for positive 
perturbation velocities and by 



for negative perturbation velocities at subsonic free stream Mach numbers. 

To check the validity of this transonic theory, comparisons have been made with 
calculations from the exact theory of Carlson using the TRANDES computer program 
(ref. 15). Mach number distributions around the blade matched well enough to justify 
use of the approximate theory, which requires two orders of magnitude less computer 
time. 
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Calculated Aerodynamic Flow Field 

Variations of streamwise perturbation velocity along the chord, as calculated for 
transonic Mach numbers from 0. 80 to 1. 04 and also for Mach numbers of 0 to 1. 15, 
are plotted in figure 6 for a 2% thickness ratio parabolic arc airfoil. As Mach number 
was increased from 0 to 0. 8, calculated perturbation velocities increased by ap- 
proximately a factor of ^1-M^j ~2 as would have been predicted by linearized sub- 
sonic theory. Further increase to 0. 90 caused a relatively larger increase of maxi- 
mum perturbation velocity so that the velocity distribution was relatively more peaked 
than that for incompressible flow. This change in shape of the velocity distribution 
agrees with that expected for the nonlinear effect of subsonic Mach number as calcu- 
lated by the Karman-Tsien compressibility correction. 

Small increases of Mach number above 0. 90 caused the development of locally 
supersonic flow followed by a shock wave. Calculated values of local velocity ratio 
upstream of midchord were approximately independent of free stream Mach number in 
the range from about 0. 91 to 0. 94. This range of Mach number was characterized by 
the change from a broad maximum of perturbation velocity near midchord to an accel- 
erating transonic flow with large perturbation velocities. 

As free stream Mach number was further increased, a general reduction in per- 
turbation velocity ratio occurred. However, the shock wave continued to move down- 
stream until it reached the trailing edge at a free stream Mach number of about 0. 97. 
The calculated distribution of perturbation velocity for a Mach number of 1. 04 resem- 
bles that for fully supersonic flow at a Mach number of 1. 15 except for having a loga- 
rithmic singularity rather than a finite value at the leading edge. This difference cor- 
responds to the presence of a detached shock wave and leading edge stagnation point at 
1. 04 but an attached shock wave at the sharp leading edge for the higher Mach number. 

These velocity distributions and the corresponding densities from eq. 37 were 
used to calculate the normalized quadrupole strength pu^/poU^ on the airfoil surface 
as shown in figure 7. In fully subsonic flow the source distribution for this airfoil is 
symmetric about midchord with sharp peaks at the leading and trailing edges and a 
broad peak near midchord. Increasing the free stream Mach number increases the 
amplitude of this broad peak. Further increase above the critical Mach number causes 
the growth of a large, strong quadrupole noise-radiating region along the aft 2/3 of the 
chord. As free stream Mach number is increased through one, the acoustic strength 
of this aft region decreases and the sharp peak near the leading edge strengthens. 
Finally, at supersonic Mach numbers large enough to maintain an attached leading edge 
shock wave, the source strength is concentrated near the leading and trailing edges. 

The integration of the quadrupole source in the vertical direction is described below; 
however, if the effective extent of the source above the surface were temporarily 
assumed to be 1 chord, then the curves in figure 7 would represent f (puVpoU^) dY. 
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This permits the shape and general level of these chordwise quadrupole distributions 
to be compared with those of the one-sided volume displacement source distribution 
H (X) /2, also plotted in the figure. For subcrltlcal Mach numbers the two distribu- 
tions are similar. This may explain observations by Schmitz and Yu (ref. 16) and 
others that the linear theories predict waveshapes like the test data but with levels 
too low. At supercritical speeds, the downstream shift of the quadrupole will result 
in a substantially different noise waveform. 

Figure 8 shows the vertical distribution of the quadrupole source for a com- 

bination of Bd» Z, and mB of interest for advanced turbopropeller applications. The 
ordinate, 


2mBBj) 



is normalized such that each curve has unit area. It can be seen that the quadrupole 
source extends farther from the airfoil at a Mach number of 0. 95 than at 0. 90 or 1. 00. 
This, coupled with the trends in pu^/poU^ on the surface as shown in figure 7, leads to 
a strong peaking of the integrated quadrupole source ilfn at transonic speeds as 
shown below. 

Figure 8 shows that the effective vertical extent of the quadrupole is about one 
chord. This means that the two-dimensional aerodynamic theory is probably adequate 
except near blade tips. It also means that for out of plane noise calculations, noncom- 
pactness in the Y direction has to be accounted for with the Y exponential in eq. 33. 

Comparison of Volume Displacement and Quadrupole Source Terms 

The aerodynamic methods described above are now used to compare the integrated 
quadrupole source with the integrated volume displacement source xj/y as func- 

tions of radius (1. e., as functions of section relative Mach number, which varies as 
radius). For the biconvex parabolic airfoil, eqs. 34 and 36 can be integrated analyt- 
ically to give 


^V(^c) 



(42) 


which is plotted in figure 9. The frequency parameter 

2mBB£) 


(43) 
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2 ^max 


which 


is a measure of noncompactness. In the limit as Bjj -^0, 


3 


b 


(multiplied by b^) is just the cross-sectional area of the airfoil. The effect of increas- 
ing the harmonic order, or increasing the chord at constant section area, is to reduce 
the noise via chordwise phase interference. The example chosen here for comparison 
with the quadrupole is the advanced turboprop mentioned above with Bd = 0. 14 and 
B = 8 . For the blade passing frequency fundamental (m = 1), figure 9 gives 

(2.24) = 0. 587 hmax/b» which is only slightly reduced due to noncompactness. 


The quadrupole source was integrated numerically in the X direction from lead- 
ing edge to trailing edge for the subsonic and transonic cases and from the bow shock 
to trailing shock for the supersonic case. Integrations in the Y direction were car- 
ried out to 20 chords, which figure 8 shows to be well within 1% of the total. The cal- 
culated ratios of ^11 to tp^y for a 2% thickness ratio airfoil section are plotted 
versus section relative Mach number in figure 10 for the 3 speed ranges. The fully 
subsonic calculations and the transonic calculations are shown to blend well in their 
range of overlap. The transonic calculations show a peak at = 0. 95 and then start 
to decay with increasing Mach number. The increase in calculated values of rpn 
starting at Mj. = 1. 05 is caused by the incorrect assumption that velocity perturba- 
tions decay along vertical lines rather than along Mach waves, so that phase cancella- 
tion is not correctly represented. The fully supersonic calculations are also shown to 
diverge at Mach numbers below about 1. 15 because the leading edge shock would not be 
attached as assumed in the theory. Thus, the transonic and supersonic calculations 
are faired together as shown by the dashed line in figure 10. 


The faired curve from figure 10 is replotted in figure 11 as 20 logj^Q 


^11 + 




which is the number of decibels added to the volume displacement thickness source by 
the quadrupole thickness source. Two features are immediately apparent. First, the 
peak value of 5. 7 dB shows clearly that quadrupole radiation is an important factor in 
high speed rotor noise. Second, the quadrupole source is impprtant only at or very 
near transonic section speeds. The acoustic results in figure 11 correspond to the 
well known aerodynamic results from thin airfoil theory in which linear sources and 
doublets give good performance predictions except near a Mach number of one. It also 
shows in a unified way why the linear source terms have been adequate for prediction 
of rotor noise at subcritical tip speeds, why linear source terms are adequate for sonic 
boom calculations (M^. > 1.5), and why linear source models fall for the transonic 
speeds of interest for advanced propellers and helicopter rotors. 


Now that quadrupole noise is recognized as a transonic flow phenomenon, one can 
ask how to suppress it. Sweep in wings is known to decrease adverse effects of tran- 
sonic flow if the effective Mach number (product of the relative Mach number and the 
cosine of the sweep angle) is less than the section critical Mach number. These favor- 
able effects of sweep have also been shown (ref. 17) to occur for rotating blades. Thus, 
by analogy figure 11 shows for 2% thickness ratio propeller sections, that the additional 
noise due to the pu^ quadrupole should be reduced to less than 1 dB by sweeping the 
blade to obtain an effective Mach number of 0. 85. 
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ROTOR NOISE CALCULATIONS 


In the preceding section, volume displacement and quadrupole sources were com- 
pared at relative Mach numbers corresponding to various radii on a blade. It remains 
to integrate these sources over the blade radius to find the net efiect of the quadrupole 
on radiated noise. This section presents noise calculations for two rotors. The first 
is a model propeller for which transonic tip speed test data are available and the 
second is a rectangular planform helicopter rotor previously studied by Farassat (ref. 
2 ). 


Transonic Propeller 

Figure 2 compared near field data from a transonic tip speed propeller with noise 
predicted using Hanson' s near field theory (ref. 4) for volume displacement and dipole 
noise surface sources. These data and predictions are replotted in figure 12 along 
with the noise predicted using a near field version of the present theory for the pu^ 
quadrupole added to the other sources. Including the qviadrupole substantially improves 
the agreement with test data, particularly in, and forward of, the plane of rotation. The 
remaining underprediction is probably caused by omission of lift effects. Lift is ex- 
pected to increase the pu^ quadrupole and will be the major contributor to the p uv 
quadrupole, which can be shown to radiate with the same directivity as the lift dipole. 
These effects are now being evaluated at United Technologies. 

Farassat Rotor 

In ref. 2 Farassat predicted thickness noise for a rectangular planform helicopter 
rotor with 10% thickness ratio parabolic arc blades and a 4% chord-to-diameter ratio. 

In ref. 4, Hanson showed that the volume displacement term from eq. 18 gave results 
equivalent to those of Farassat. Figure 13 shows the radial distribution of the addi- 
tional source strength due to the pu^ quadrupole for a tip rotational Mach number of 
1.1. As in figure 11, the ordinate is 20 logj^Q |(^11 + thick 

blade, the peak occurs at a lower section speed (Mj. = 0. 88) as would be expected be- 
cause of the reduced critical Mach number. However, the peak increase of 6. 2 dB is 
nearly the same as for the thinner airfoil (as well as intermediate thicknesses). Since 
this result was a surprise, it was verified by using an exact inviscid-flow transonic 
airfoil computer program (ref. 15). It was found that, while the surface values of the 
quadrupole grow roughly as (huiax/b)3/2 as predicted by transonic similarity laws, 

the decay rate with vertical distance increases roughly as (hjjjg^x/^) ^ • The net effect 
is that the peak value of the quadrupole source grows approximately as (h ^y,y A»), 

like the surface source. However, as can be seen by comparing figures 11 and 13, the 
quadrupole is significant compared to the volume displacement over a wider range of 
relative Mach niunbers for the tMcker airfoil. 

The effect of the quadrupole integrated over the full radius of the Farassat blade 
has been checked with eq. 32. For the first several harmonics the noise increase is 
3. 8 to 4. 0 dB in the plane of rotation. However, many more harmonics would have to 
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be evaluated to find the overall effect on the waveform because of the impulsive nature 
of the sound from this blade. 


CONCLUSIONS 

An analysis has been presented for the far field harmonic thickness noise of a 
nonlifting rotor. The acoustic radiation equations are essentially exact except for 
the thin blade approximation. The pu^ quadrupole has been evaluated using two- 
dimensional aerodynamic theory valid throu^ the transonic speed range. Compari- 
son of this quadrupole with the volume displacement sources has led to the following 
conclusions. 

1. The pu2 quadrupole is an important thickness noise source at trsinsonic blade 
section speeds. For any thickness ratio, the maximum contribution is roughly 
6 dB above the volume displacement thickness noise, 

2. For thin propeller blades the quadrupole becomes negligible outside the transonic 
speed range, which explains why linear source acoustic theories are successful 
for subcritical propellers and for sonic boom calculations. 

3. As thickness ratio is increased, the quadrupole contributes significantly in the 
range of transonic section speeds found in current helicopter designs. 

4. The comparison of theory with test data for a transonic tip speed propeller was 
substantially improved by adding the pu^ quadrupole to the linear source terms. 

5. Quadrupole radiation should be reduced to the negligible point if the blades are 
swept so that their effective Mach number is subcritical. 
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APPENDIX 


Blade-Fixed Coordinates for Time Domain Calculations 

As was explained in ref. 4, the time domain equations can be transformed to 
blade-fixed coordinates by substituting 

Tj = r + Ut - MyR (44) 

For example, the pu^ quadrupole pressure from eq. 17 can be written 

Pu<^> - 4 k -ft / ff Til (y + Ut - M,E> d| dy dr. ( 45 ) 

Differentiation of eq. 44 gives 


0R 

drj = dy - Mp — dy 
^ ay 


(46) 


8R 


But — = Mor is the Mach number of the source relative to the observer so that 

■t rv'v 


dy can be expressed dy = drj/ jl -Mor j . Substitution into eq. 45 gives 

1 d f 1 


1-Mr. 


d (vol) 


(47) 


where d (vol) = d^ dr) dro . Equation 47 would be practical for numerical computation 
for subcritical rotors where 1 - M^j. never goes to zero and there are no shock waves 
to cause discontinuities in pu^ . 
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Figure 1.- Transonic tip speed propeller tested in United Technologies Research 
Center Acoustic Wind Tunnel. Model diameter = 0.62 m (2.04 ft). 
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Figure 6.- Calculated chordwise velocity distributions for 2% thickness ratio 

parabolic arc airfoil. 
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QUADRUPOLE SOURCE TERM ON AIRFOIL SURFACE, P^^IP 



SECTION RELATIVE MACH NUMBER, M^ 


lORDWISE DISTANCE, X 


;e strength along the surface of a 2% 
airfoil at transonic speeds. 
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NORMALIZED CHORDWISE INTEGRAL 



Figure 8.- Vertical distribution of quadrupole source strength. 90% points 
indicate value of Y which includes 90% of the quadrupole integral. 



Figure 



2m B Bq 

NONCOMPACTNESS PARAMETER, O)- = 

Z 

9.- Parabolic arc thickness distribution H(X) and its frequency domain 

counterpart • 
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NOISE ADDED BY QUADRUPOLE -dB 
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Figure 11.- Increase of blade thickness sound pressure level caused by 

including quadrupole noise. 




Figure 12.- Increase in noi 

data. 
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